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ABSTRACT 

Given its extraordinary spatial resolution and sensitivity, the projected Next Generation Space Telescope (NGST) 
is likely to detect a large number of high-redshift QSOs lensed by spiral galaxies. Using realistic models for the 
QSO and spiral populations, we calculate the expected number density of detectable QSOs multiply imaged by 
spirals, and investigate the influence of various evolution effects on that number density. It is shown that NGST 
will probably find of order ten lensed QSOs per square degree at 26th magnitude in the V and L bands, and that 
various observable quantities like the total number density of lensed QSOs in these two bands, the ratio between 
the number densities of lensed QSOs in the V and L bands, the fraction of QSOs with more than two images and 
so forth can be used to constrain the evolution of the QSOs, the spirals, the dust in spirals, and the masses of spiral 
disks. 
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While strong lensing of QSOs by elliptical galaxies is a well- 
evolved subject (e.g. Turner, Ostriker & Gott 1984; Fukugita 
& Turner 1991; Maoz & Rix 1993; Kochanek 1996a,b; Kee- 
ton, Kochanek & Seljak 1997; Falco, Kochanek & Munoz 1998; 
Keeton, Kochanek & Falco 1998), interest in lensing by spiral 
galaxies has only started to develop over the past few years 
rMaller, Flores & Primack 1997; Keeton & Kochanek 1998; 
^ Bartelmann & Loeb 1998; Blain, M0ner & MaUer 1999). This 
Q^was legitimate for the interpretation of QSO surveys undertaken 
JL at angular resolutions of order 1". Once angular resolutions 
^ on the order of 0.1" can be achieved, a substantial portion of 
the lens population may be spiral galaxies. If observations at 
high angular resolution can be combined with high sensitivity, 
• • numerous QSOs multiply imaged by spiral galaxies can be ex- 
. ^ pected to be observed, and such observations can provide de- 
k>( tailed information on the evolution of spiral galaxies and their 
''-^ properties out to fairly high redshift. 

This paper addresses the question of how many QSOs lensed 
by spiral galaxies can be expected to be observed by telescopes 
like NGST, whose angular resolution exceeds 0.1", and which 
reach a sensitivity of order one nano-Jansky in the visible and 
near-infrared wave bands. Going to such observational limits, 
QSOs lensed by spiral galaxies should be discovered in substan- 
tial numbers. Many of the multiple-image systems produced by 
spiral galaxies will straddle near-edge-on galactic disks. The 
statistics of image separations and time delays will therefore al- 
low to determine properties of spiral disks at moderate and high 
redshifts. The wavelength regime between the visible and the 
near infrared is particularly interesting for such investigations, 
since the importance of dust changes dramatically across this 
wavelength range. While multiple QSO images, in particular 
near galactic disks, are substantially obscured by dust in the 
V band, they are almost unaffected in the L band. Comparing 
observations taken in these two wavebands, it will therefore be 
possible to quantify the amount and the distribution of dust in 
distant spiral disks. 

The very faint flux limits envisaged will bring high-redshift 
QSOs into view. The QSO population at high redshifts is un- 
certain. Observations show a quick drop of the QSO number 
density beyond redshifts of 4 — 5 (Pei 1995). This may reflect 
rapid QSO evolution at these redshifts, but also incompleteness 
of the QSO samples observed. We therefore investigate two al- 



ternatives, assuming either that the observed QSO population is 
complete, or that it is to be extrapolated towards high redshifts. 

It is also important for this study to reliably estimate QSO K- 
corrections in the visible and the near infrared. For this purpose, 
we model the dominant features of QSO spectra together with 
the Lyman-a forest, and determine /T-corrections by convolving 
the redshifted model spectrum with the appropriate filter func- 
tions. 

Finally, the contribution of disks to lensing by spiral galax- 
ies is crucial for two reasons. First, projected disks at large 
inclination angles reach high surface mass densities and affect 
multiple-imaging cross sections substantially. Second, dust in 
these disks demagnifies images and thereby alters the statistics 
of observed image properties. 

We proceed in the following way. Section 2 describes the 
computation of QSO /^-corrections and summarises the model 
for the QSO luminosity function. Our model for the spiral galax- 
ies is introduced in Section 3. Section 4 briefly describes the 
numerical method. Section 5 presents the results, and we finish 
with the conclusions in Section 6. 

Throughout, we adopt a CDM cosmological model with Q. = 
0.3, £2a = 0, and a Hubble constant of Hq = 70kms"' Mpc"'. 
The CDM spectrum is normalised such that the local abundance 
of rich galaxy clusters is reproduced (Eke, Cole & Frenk 1996; 
Viana & Liddle 1996). 

2. QSOs 

2.1. QSO spectra and the K-correction 

A reliable estimate for the QSO /iT-correction is very important 
for the present study. We therefore have to construct a synthetic 
QSO spectrum representing the prominent features of a typical 
QSO spectrum. This spectrum consists of the underlying QSO 
continuum, the hydrogen emission lines of the Lyman, Balmer, 
and Paschen series, the corresponding Balmer and Paschen con- 
tinua, the two-photon continuum, the 3000 A bump mainly con- 
tributed by iron (Fe II) emission lines, and emission lines of sev- 
eral other elements like helium, carbon, nitrogen, oxygen and so 
forth. 

We adopt the QSO continuum emission from Elvis et 
al. (1994), who tabulate the median QSO continuum in the fre- 
quency range 10^ Hz < v < lO'^Hz. The spectrum on top of 
the broad continuum is modelled as in Grandi (1981, 1982). It 



contains the Balmer- and Paschen series up to order 28 with 
line strengths adopted from Brocklehurst (1971), exponential 
Balmer- and Paschen continua, the two-photon continuum as 
given by Spitzer & Greenstein (1951), 79 emission lines of 
Fe II to reproduce the 3000 A bump, and several other emis- 
sion lines of helium, carbon, nitrogen, oxygen, magnesium, and 
silicon. The electron temperature is assumed to be 12, 500 A', 
the electron density is 10^ cm^^, and the lines are broadened 
by Lorentzian profiles with an FWHM velocity of 4, 500km s^ ' . 
Some additional lines tabulated by Netzer & Davidson (1971) 
are added. Following Grandi (1982), the emission lines are nor- 
malised such that the equivalent width of the Hp line against the 
continuum is 100 A. Finally, the number density and equivalent 
widths of Lyman-a forest lines are modelled following Murdoch 
et al. (1986). A synthetic spectrum of a OSO at redshift z = 4 
obtained in this way is displayed in Fig. (P. 
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Fig. 1 . — Synthetic QSO spectra at redshifts z = and z = 4 (top and 
bottom panel, respectively), including strong Lyman-a forest lines. 
Wavelengths range from 900A to 36, 000 A. The ordinate is scaled to 
arbitrary units. The most prominent line in both spectra is Lyman-a. 

The synthetic QSO spectrum is then convolved with the re- 
sponse functions of the Johnson V and L bands to obtain the K- 
correction. The peak transmission wavelengths of the response 
curves are at 0.555 /jm for the V band and at 3.6/jm for the L 
band. The /iT-corrections in these bands are plotted in Fig. (^. 

The A'-correction in the V band is negative from redshift zero 
to just above redshift 3, where it rises steeply. Negative K- 
corrections are due to bright emission lines moving into the V- 
band response function, while the steep increase is caused by 
Lyman-a forest absorption. In the L band, the /iT-correction 
ranges between ±1 magnitudes. The shallow peak in the L- 
band /iT-correction is due to the shallow valley in the QSO spec- 
trum beyond w 8, 000 A. The pronounced minimum atz~ 4.4 is 
caused by the Balmer-a line's coming into the L-band response 
function. 
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Fig. 2. — AT-corrections as functions of QSO redshifts in the V and L 
bands (solid and dashed lines, respectively), computed from synthetic 
QSO spectra convolved with the V and L response functions. 



2.2. QSO luminosity function at high redshift 

Haiman & Loeb (1998) recently provided a model for the QSO 
luminosity function at high redshifts. The model rests on four 
simple and intuitive assumptions: (i) black holes form in all 
dark-matter haloes above a mass limit Mo; (ii) the black-hole 
mass Mbh is a certain constant small fraction e of the halo mass; 
(iii) black holes produce QSOs shining at the Eddington lumi- 
nosity Ledd = i'edd,© at birth, and the QSO luminosity is 
on average a monotonically decreasing function Mbh /(O of time 
t after birth; and (iv) the mass function of dark-matter haloes is 
given by the Press-Schechter (1974) mass function N{M). 
The QSO luminosity function is then 

4>(L,z) = fdz' r dM^^5[L-Mbh/(f)] • (D 

The derivative of the Press-Schechter mass function with respect 
to z gives the birth rate of haloes and thus the birth rate of black 
holes, and the delta function selects those QSOs shining with 
luminosity L at time t after their birth. 

The redshift derivative of the Press-Schechter mass function 
contains a negative contribution because of halo merging. It 
therefore underestimates the halo formation rate. However, the 
merger probability at high redshift is low, so that the underesti- 
mate is unimportant (see also Sasaki 1994). 

The black-hole mass fraction e is kept constant here for sim- 
plicity. Introducing a moderate scatter in e (cf. Magorrian et 
al. 1998; Van der Marel 1999) affects the black-hole mass func- 
tion only weakly (Haiman & Loeb 1999). 

Haiman & Loeb used an exponential function for /(f) and 
assumed Mq = 10^ Mg [(1 +z)/10]-^/^ for the minimum halo 
mass. The model has then only two free parameters, the black- 
hole mass fraction e and the time constant to of /(f). Haiman 
& Loeb showed that e « 5.4 x 10"^ and fo « 6.3 x 10^ yr pro- 
vide an excellent fit to the observed QSO luminosity function at 
intermediate redshifts (Pel 1995). 



Carrying out the redshift integral in (M leads to 
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Here, z* is the redshift at which a QSO in a halo of mass M 
needs to be born in order to reach luminosity L at redshift z, 
and is the time elapsed between z* and z- Note that another 
mass limit on the haloes is implicit here. The QSO cannot shine 
with higher than Eddington luminosity. The halo mass therefore 
needs to exceed both Mq and (Mq/e) (L/Ledd,©)- 

We repeat the derivation of Haiman & Loeb's luminosity 
function here because numerical evaluation reveals that suffi- 
cient approximations can be made to render the remaining in- 
tegral in analytically tractable. This is convenient because 
the accurate numerical evaluation of ^ requires substantial care 
because of the large dynamic range of the Press-Schechter mass 
function. 

We introduce two approximations. The first is that is typi- 
cally very short compared to the Hubble time, so that w z is 
a very good approximation. The second is that the power spec- 
trum P{k) of the dark-matter density fluctuations can be written 
as a power law, 

P{k) oc k" . (3) 

This approximation is well satisfied in the range of scales k 
corresponding to the relevant dark-matter halo masses, where 
« w —2.4. The first approximation removes the necessity to nu- 
merically evaluate z*. The second approximation allows to write 
the Press-Schechter mass function as 
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where p is the mean cosmic density today, is the nonlinear 
mass scale today, m = M/M^, is the halo mass in units of M*, 
D{a) is the growth factor of density perturbations as a function 
of the scale factor a, and a is related to n through 



a ; 



« + 3 



(5) 



With (0), the remaining mass integral in can be carried out. 
The finalresult for the luminosity function is 



X [2r(i + p,rt))-r(p,A<o)] , (6) 

where we introduced the abbreviations 

a- 1 



and 



2a 



..2a 



As mentioned before, /io is determined through by 

i'edd.oE 



Mq — max Mo, 



(7) 



(8) 



(9) 



and m = M/M^, . For low QSO luminosities, the term in brackets 
in (||) is constant in luminosity, and <I>(L.z) °^ L^^. 



For n = —2.4, a = 0. 1 and P = —4.5. In an Einstein-de Sitter 
universe, D{a) = a and dz/df = Hoa^^l^, and thus 
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(10) 



Equation ( |lO| ) shows that the QSO luminosity function <I>(L,z) 
is a very steep function of redshift. Ignoring the redshift de- 
pendence of Mo, <I>(L,z) scales with a approximately as a^23/2^ 
strongly depending on the local slope of the dark-matter power 
spectrum. 

2.3. QSO number counts 

Haiman & Loeb's model is valid at high redshifts. At low 
redshifts, the QSO luminosity function can be assumed to be 
sufficiently well observed. In order to test the impact of the 
theoretical extrapolation on the expected counts of strongly 
lensed QSOs, we use two alternative QSO luminosity func- 
tions, namely either Pei's (1995) fit to the observed QSO number 
counts alone, or Pei's function augmented by Haiman & Loeb's 
(1998) theoretical model. In the latter case, we construct the lu- 
minosity function <P{L,z) from the observed function out to red- 
shift 3, from the theoretical model beyond redshift 4, and from 
a redshift-dependent linear combination of the two between red- 
shifts 3 and 4. For simplicity, we shall refer to these two lumi- 
nosity functions as the "observed" and the "theoretically extrap- 
olated" one below. 

Using the /T-correction derived in § 2.1, the QSO luminosity 
function can now be converted to number counts per flux and 
redshift intervals dS and dz. 



d^NjS.z) 
dSdz 



dSdz. 



(11) 



3. Spiral GALAXIES 



3.1. Evolution of spirals 



We assume that the spiral population has a Schechter (1976) lu- 
minosity function, i.e. the number density n follows 



dn{l) 
dl 



- n* r exp(— /) 



(12) 



where / — L/L^, is the luminosity in units of L^,, n^, — 1.5 x 
10^^/i^Mpc^^ is the present-day number density, and v = 
—0.81 (Marzke et al. 1994). We further assume that the cir- 
cular velocity scales with I according to the Tully -Fisher (1977) 
relation. 



— = / 

V, 



(13) 



with a = 4. We adopt two scenarios for the redshift evolution of 
the spiral population, namely either constant comoving number 
density, or a number density evolving according to the Press- 
Schechter mass function. 

Mo, Mao & White (1998) suggested a model for the evolution 
of spiral disks resting upon four simple and intuitive assump- 
tions. These are: (i,ii) The disk mass and angular momentum 
are fixed fractions of the spiral mass and angular momentum, 
(iii) The radial disk profile is exponential, and the disk is cen- 
trifugally supported, (iv) The disk is dynamically stable. 



These assumptions yield a set of simple scaling relations for 
the disk properties. Disk mass Md(z) and radius scale with 
the inverse Hubble parameter H{z), 



H{z) 



Ho 



-1 



(14) 



and consequently the disk surface mass density Ed{z) scales with 
the Hubble parameter. 



Ho 



(15) 



We use either constant disk parameters below or evolve them 
according to eqs. (|l3) and (ll). 

Finally, we assume that the disk radius scales with lumi- 
nosity as 

rd-/'/^ (16) 

3.2. The lens model 

Our mass model for the spirals consists of two components, a 
halo and a disk. They are modelled as described in Keeton & 
Kochanek (1998) and Bartelmann & Loeb (1998). We refer 
the reader to these papers for details. In essence, both com- 
ponents are described as isothermal, arbitrarily flattened, ellip- 
soidal mass distributions, characterised by a core radius s and 
an oblateness q, defined as the ratio between the short and the 
long axes of an oblate ellipsoid. We can then compose the lens- 
ing potential i|/ of a spiral with a maximal disk of three building 
blocks. 



¥(«d,?d) - v('-d,?d) + \|/(ih, 1) 



(17) 



The first two terms are the lensing potential of the disk, truncated 
at radius rd, the third is the lensing potential of the halo. For sub- 
maximal disks with a disk fraction < /d < 1, this is changed 
to 

¥ = /d¥max + (l-/d)¥K,l) , (18) 

i.e. the truncated maximal-disk model is scaled down by /d and 
embedded in a spherical halo of unchanged mass. We therefore 
have five parameters in total for the mass model, namely the core 
radius of the disk id, the core radius of the halo s^, the oblateness 
of the disk q^, the disk fraction /d, and the asymptotic circular 
velocity Vc. 

For an galaxy, we assume a circular velocity of Vc* = 
220km s^' and a disk radius of rd* — 8/i~^kpc. They scale with 
spiral luminosity as given by eqs. (Oh and (n6h. 

Flat rotation curves require a halocore raHius of w 0.72 rd 
(Keeton & Kochanek 1998), and for the disk core radius we as- 
sume id = 0.2/i^'kpc. The disk fraction is either set to 1 or 
0.5. 

Up to a constant, the light travel time to the observer from a 

source at angular position P through an image position 9 is given 
by 
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(19) 



where z' and z are the lens and source redshifts, respectively; 
D{z' ,z) is the angular-diameter distance between redshifts z' and 

z, and D{z) = D{Q,z). The lens equation relates and P through 



p = e-a(e), 



(20) 



where a(9) = V\)/(9) is the deflection angle. The time delay 
between images is therefore determined by the lensing potential, 
its gradient, and the geometry of the lens system. 



3.3. Gas and dust 

The disk is further assumed to encompass a double-exponential 
neutral-hydrogen disk with a scale radius given by rd, a scale 
height of ^d'"d, and a central column density of 11.25 x 
10^" cm"^ (Broeils & van Woerden 1994). Dust follows the 
neutral-hydrogen distribution with a present-day dust-to-gas ra- 
tio of 0.01 (Whittet 1992). The dust is assumed to be a mixture 
of silicates and graphites, for which the extinction curve as a 
function of wavelength was derived by Draine & Lee (1984). 
The dust-to-gas ratio is either assumed to be constant with red- 
shift, or to evolve in proportion to (1 +z)^^ (Pei, Fall & Bech- 
told 1991). 



3.4. Summary of the spiral model 

Having fixed the disk radius, the disk core radius, the halo core 
radius, and the asymptotic circular velocity, and assuming the 
scalings of the circular velocity and the disk radius with lumi- 
nosity, the mass model for the spirals has two more free param- 
eters, namely the disk fraction /d and the disk oblateness q^. 

The number density of spirals is then either assumed to be 
constant or to evolve according to the redshift evolution of the 
Press-Schechter mass function in the appropriate mass range. 
The parameters of individual spirals are either assumed to be 
constant with redshift, or to evolve according to the model by 
Mo et al. (1998). The dust-to-gas ratio in the spiral disks finally 
is either assumed to be constant or to evolve proportional to ( 1 + 

We focus on a fiducial model with disk oblateness q^ — 0.2, 
disk fraction /d = 1, and evolving disk parameters, number den- 
sity, and dust. 

3.5. Visibility of faint QSO images 

Many QSO images lensed by spiral galaxies will appear close 
to spiral disks. In order to be detectable, they need to be suf- 
ficiently brighter than the disk itself. We can estimate a maxi- 
mum QSO magnitude in the following way. Freeman's (1970) 
law asserts that the central surface brightness of nearby spi- 
rals is w 21.5magarcsec^^ in the B band, corresponding to 
/o,v ~ 20.8 mag arcsec^^ in the V band. Ignoring spectral evolu- 
tion, this changes to 



/o,v(z) =/o,v-2.5(a-3)logio(l+z) , 



(21) 



at higher redshift, where a « — 1 is the power-law exponent of 
the spectral-energy distribution. The most probable lens redshift 
for high-redshift QSOs is z w 0.5 — 0.6, where the central surface 
brightness is /o,v(0.6) ~ 22.6 mag arcsec^^. 

Typical image separations from the lens centre are of order 
the Einstein radius, or roughly 1 h^ ' kpc. Assuming a disk scale 
length of « 8/;^' kpc and an exponential light profile, the disk 
surface brightness near image locations has dropped by a factor 
of w 0.9 to /v(0.6) w 22.7 mag arcsec-2 nqsT is planned to 
have a pixel size of order 5 x lO^^arcsec^, hence the per-pixel 
disk magnitude is of order 3 1 mag. Assuming that a QSO im- 
age covers a few pixels only, it seems realistic to expect that 
QSO images brighter than « 28 — 29 mag can safely be detected 
on top of a spiral disk at redshift « 0.5 — 0.6. However, many 
multiple QSO images will straddle near-edge-on disks, whose 
surface brightness drops rapidly away from their symmetry line. 

This estimate applies to the V band. Spiral disks are likely 
to be brighter in the near infrared, where QSOs are intrinsi- 
cally fainter. However, QSO host galaxies are substantially 



more prominent in the near infrared and give rise to character- 
istic lensing features like arcs and rings (Kochanek et al. 1999), 
which help in detecting lensed QSOs against spiral disks, and in 
distinguishing them from surface-brightness fluctuations. Con- 
fusion of lensed QSO images with bright spots of stellar ap- 
pearance like OB associations or other surface-brightness fluc- 
tuations can also be disentangled through colour information or 
highly resolved X-ray observations. The order of magnitude of 
the previous estimate should therefore also hold in the L band. 

4. Method of computation 

We wish to compute the number of QSOs above a certain flux 
limit S lensed by spiral galaxies such that the images have certain 
properties Q. Let 



dSdz 



-dSdz 



(22) 



be the number of unlensed quasars with flux within dS of S and 
redshift within dz of z- Let further P(/j,z;2) be the probability 
for a quasar at redshift z to be magnified by a factor > /j and 
lensed with image properties Q. The number of lensed quasars 
in the redshift interval zi <z<Z2 above the flux limit S is then 



N{S) 



dS' rdzp(^,z;Q 



d^N{S,z) 
dSdz 



iS',z). (23) 



The lensing probability P{iJ,z',Q) is conveniently computed 
in the cross-section approach. The cross section a{iJ,z' ,z,l',Q) 
is defined as the area on the source sphere at redshift z in which 
a source must fall in order to be lensed with image properties Q 
and magnified by > /J by an individual galaxy with luminosity 
/L* at redshift z'. The lensing probability is then 



X J^o(A<,z',z,/;e), (24) 
dz 



where (1 + z!) n(l.,z!) is the number density of galaxies with 
luminosity I at redshift z', [dy(z')/dz]dz' is the proper volume 
of a spherical shell of width dz' and radius z', and D{z) is 
the (angular-diameter) distance from the observer to the source 
plane. 

It is important to note that total image magnifications are not 
only determined by lensing in our model, but also by dust extinc- 
tion in the galactic disks. Total magnifications therefore range 
within < /J < oo rather than 1 < /j < oo, and therefore S' in ( P3| ) 
ranges within Q < S' < °°. 

The scaling relations (13) and dl^ ) assure that the cross sec- 
tion o scales with the squared Einstein radius. The Einstein ra- 
dius in turn scales with the squared circular velocity, and hence 
with Assuming the scaling (|l6|) of the disk radius, the disk 
parameters scale with galaxy lummosity / in the same way as 
the Einstein radius. This allows to immediately perform the l- 
integral in (|24|). The remaining redshift integral needs to be per- 
formed numerically. 

We are left with the determination of the cross sections o. 
Since the galactic disk defines a symmetry plane randomly in- 
clined with respect to the line-of-sight, this involves averaging 
the cross sections over inclination angles. For fixed inclination 
angles, cross sections are computed by placing sources on a grid 
in the source plane, finding all images of each source, computing 
their magnifications and classifying the image properties, and 
finding the area in the source plane in which sources are magni- 
fied by > /J and imaged with properties Q. 



Technically, we increase numerical efficiency by placing the 
sources on an adaptive grid in the source plane in order to 
increase numerical resolution where needed, i.e. near caustic 
curves. Rough image positions are then computed in a first 
coarse scan in the lens plane, followed by a fine scan to re- 
fine image positions and to determine accurate image properties. 
This technique of having adaptive grids in both the source and 
the lens planes increases computational speed by more than an 
order of magnitude compared to the straightforward fixed-grid 
approach. Different levels of grid refinement in the source plane 
are statistically accounted for by assigning appropriate statistical 
weights to the sources. 

Having found the cross sections o, we can evaluate the lens- 
ing probability (Q). With the quasar number counts taken from 
§ 2.2, the number of quasars above flux limit S and imaged with 
properties Q is then computed evaluating (^3]). 

5. Results 

Multiply imaged QSOs can be recognised as such if the angu- 
lar separation A9 of the images exceeds a certain threshold, and 
if the flux ratio r between the images does not exceed a cer- 
tain dynamic range. We adopt A9 > 0.1" and r < 20 throughout 
this section; i.e., with "lensed QSOs" we mean multiply imaged 
QSOs whose images satisfy these criteria. In case of three or 
more images, we take A9 and r between the two brightest im- 
ages. 

The number densiU' of lensed QSOs per square degree on the 
sky is shown in Fig. ^as a function of total apparent QSO bright- 
ness. The figure has three panels, and each panel contains four 
curves. Panel (a) shows results for the fiducial model. It assumes 
a spiral population with evolving number density and evolving 
disks, containing evolving dust. Dust evolution is neglected in 
panel (b), and the evolution of the spiral population is ignored in 
panel (c). 

The solid and dotted curves show the number density based 
on the observed QSO luminosity function in the V and L bands, 
respectively. The L-band (dotted) curve falls below the V-band 
curve by approximately one order of magnitude because QSOs 
are intrinsically fainter in red than in blue light. The remaining 
two curves (short- and long-dashed) are again V- and L-band 
results, now based on the theoretically extrapolated QSO lumi- 
nosity function by Haiman & Loeb. The changed QSO lumi- 
nosity function affects the V-band counts (solid vs. short-dashed 
curves) only by a very small amount, while it has a very signifi- 
cant effect on the L-band results. The reason is the K-correction. 
Although the theoretically extrapolated QSO luminosity func- 
tion extends to much higher redshifts than the observed one, the 
high-redshift QSOs are almost entirely suppressed in the V band 
because of the strong K-correction setting in at redshifts > 3.4. 
In the L-band, however, the K-correction remains small or neg- 
ative even at high redshifts (cf. Fig. 0). Consequently, the huge 
reservoir of high-redshift QSOs predicted by the theoretically 
modelled QSO luminosity function is not at all suppressed in 
the L-band, where the different assumptions for the QSO lumi- 
nosity function result in number densities differing by approxi- 
mately one order of magnitude at apparent magnitudes of w 28. 
Faintward of about 26.5 th magnitude, the L-band counts exceed 
the V-band counts. 

The neglect of dust evolution in panel (b) has a strong effect 
on the V-band counts, and a substantially smaller effect in the L 
band. V-band counts are reduced by factors of 3 — 4, regard- 
less of the assumed QSO luminosity function. L-band counts 
based on the observed QSO luminosity function are lowered by 
w 25% compared to evolving dust, and approximately halved 
if the theoretically extrapolated QSO luminosity function is as- 
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Fig. 3. — Number of quasars on the sky per square degree, lensed by spiral galaxies such that there are two or more images with image separation 
> 0.1" and flux ratio < 20. Each panel has four curves. Solid and short-dashed curves: V band, QSO luminosity function fit to the observed 
number counts only, or extrapolated by the theoretical model; dotted and long-dashed curve: likewise for the L band. The three panels show: (a) 
spiral population and dust evolve; (b) only spiral population evolves; (c) only dust evolves. 



sumed. The reason is of course that with increasing QSO red- 
shift, the redshift of the most efficient lenses also increases, so 
that QSO light in the observer's L band is more strongly affected 
by dust in the lenses. This leads to a relative enhancement of the 
L-band compared to the V-band counts. 

When the spiral number density is assumed to be constant 
in comoving space, as in panel (c), the lensed QSO counts are 
reduced by almost an order of magnitude in both bands. This 
is due to the fact that spirals merge to form ellipticals towards 
redshift zero, hence the spiral number density is expected to in- 
crease towards moderate redshifts. If this increase is neglected, 
there are fewer lenses available and therefore fewer lensed QSOs 
observable. 

Variation of other parameters has very little effect on the 
counts of lensed QSOs. Lowering the disk mass (by reducing /d 
from unity to 0.5), or neglecting the expected change in spiral 
properties with redshift, does not change the lensed QSO counts 
significantly. 
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Fig. 4. — Time-delay distribution for images split by > 0.1" with flux 
ratio < 20 and total apparent magnitude < 25 in the V and L bands. 
Three curves are shown per panel. Solid curve: fiducial model; dotted 
curve: no scale evolution; dashed curve: disk fraction reduced to = 
0.5. 

Figure Q shows the distribution of time delays between the 



two brightest images of lensed QSOs. The figure has two panels 
for the two filter bands, as indicated. Each panel shows three 
curves. The solid curve is the distribution for the fiducial model, 
in which spiral evolution is taken into account, and the disks are 
maximal. The curve peaks at time delays of order 0.02 days, 
or approximately half an hour It goes through a local mini- 
mum at w 0.2 days, and rises again for a broad and shallow 
peak centred on w 1 day. This bimodal distribution reflects the 
two components of the mass distribution. Time delays of order 
one day occur between images whose splitting is dominated by 
the halo, while images straddling a near-edge-on disk have sub- 
stantially smaller time delays. This is further demonstrated by 
changing the disk contribution to the lens model. When disk 
evolution is neglected, as for the dotted curve, the time-delay 
distribution peaks more sharply at small time delays, and the 
secondary maximum indicative of the halo disappears. Con- 
versely, the secondary maximum is emphasised relative to the 
maximum at small time delays when the disk fraction is reduced 
to /(J = 0.5, as for the dashed curve. 

The secondary maximum in the time-delay distribution is 
more pronounced in the V than in the L band. The reason is 
the dust in the galactic disks. Images straddling near-edge-on 
disks, which are responsible for the peak at low time delays, are 
strongly suppressed by dust in the V band, but easily seen in the 
L band. The primary maximum in the time-delay distribution is 
therefore more pronounced in the L than in the V band, so that 
the secondary maximum is relatively lowered. 

The results in Fig. Q are based on the theoretically extrapo- 
lated QSO luminosity function, but this assumption is insignifi- 
cant in this context. Adopting the observed luminosity function 
leads to virtually identical results. 

Figure ^ shows the fraction of multiple-image systems with 
three or more visible images. As before, images are consid- 
ered visible if they are separated by > 0.1" and span a dynamic 
range of < 20. Although these systems include such with four 
or more images, we refer to the fraction as "triple fraction" for 
simplicity. The figure also has two panels for the two filter bands 
considered. 

Both panels have four curves. The solid curve represents the 
fiducial model, which assumes evolution of the spirals and of 
the dust in the spiral disks. There are w 30% image systems 
with more than two images in the V band, and « 40% in the 
L band. It is again straightforward that there should be more 
such systems in the L rather than in the V band. Dust in the 




Fig. 5. — Fraction of multiple image systems with three or more ob- 
servable images for the V and L bands. Solid curve: fiducial model; 
dotted curve: no dust evolution; short-dashed curve: no scale evolu- 
tion; long-dashed curve: disk fraction reduced to = 0.5. 



Fig. 6. — Number d/V /dz of lensed QSOs with total apparent magnitude 
< 25 per square degree and unit redshift. Solid curve: V band, observed 
QSO luminosity function; dotted curve: L band, observed QSO lumi- 
nosity function; short-dashed and long-dashed curves: theoretically ex- 
trapolated QSO luminosity function, V and L band, respectively. 



disks suppresses images and therefore converts part of the triple 
systems to doubles, and this effect is much more important in 
the V than in the L band. 

Neglecting dust evolution leads to the dotted curves. The 
triple fraction is reduced to just below 10% in the V band and to 
« 30% in the L band. Obviously, the reduction compared to the 
model with evolving dust is more severe in the V than in the L 
band. 

The long-dashed curves are obtained after reducing the disk 
fraction from unity to = 0.5. Generally, the triple fraction is 
reduced, but now the reduction is more severe in the L than in 
the V band. The reason for this is that dust suppresses triple 
systems in the V band much more strongly than in the L band. 
Therefore, the reservoir of triple images in the L band is much 
larger, so the relative effect of a lower disk fraction is stronger 
in the L band. 

The strongest effect, however, is achieved by neglecting the 
evolution of spiral disks, as shown by the short-dashed curve 
in Fig. I Neglecting this aspect of the evolution leads to more 
extended and dominant disks at moderate redshifts. Such disks 
cause pronounced caustic curves with narrow, extended spikes 
following the projected disks. These extended caustics increase 
the multiple-image cross section of spirals considerably. As 
Fig. H shows, the triple fraction is increased to just below 60% 
in the V band, and to almost 85% in the L band. 

The results shown in Fig. ^ are based on the theoretically ex- 
trapolated QSO luminosity function. Generally, the V-band re- 
sults are unaffected by the choice of luminosity function, and the 
L-band results are generally slightly higher if the observed QSO 
luminosity function is used instead. The reason is highlighted 
by Fig. I 

Figure |6| shows the redshift distribution of lensed QSOs in 
the V and L bands which are brighter than 25 th apparent mag- 
nitude. The solid and the dotted curves show the V-band dis- 
tribution with the observed and the theoretically extrapolated 
QSO luminosity function, respectively. The two distributions 



are very similar, and they peak just above z = 3. The distribu- 
tion in the L band based on the observed QSO luminosity func- 
tion (long-dashed curve) shows essentially the same behaviour, 
except that its amplitude is lower because of the lower L-band 
QSO counts at 25 th magnitude. Assuming the theoretically ex- 
trapolated QSO luminisity function, however, broadens the peak 
and shifts it to z « 5 — 6. The substantially higher QSO red- 
shifts in the L band lead to more important dust extinction, and 
therefore to lower triple fractions, if the theoretically extrapo- 
lated QSO luminosity function is assumed. 

In addition, this explains the different trends with apparent 
QSO magnitude of the curves in Fig. ^. While the triple fraction 
is essentially constant across the plotted range of magnitudes in 
the V band, it decreases for fainter QSOs in the L band. This 
is because the faint QSO population originates from a broader 
redshift range in the L than in the V band, and the importance of 
dust extinction increases with QSO redshift. 



6. Conclusions 

We modelled strong lensing of high-redshift QSOs by spiral 
galaxies, aiming at future surveys with high angular resolution 
and high sensitivity in a wavelength range between the visible 
and the near infrared. The main question was, how many mul- 
tiply imaged QSOs are expected to be seen, and how do their 
number and their image configurations depend on the properties 
of the lensing spiral population? 

Our mass model for the spiral galaxies consists of a spherical 
halo and a disk, which also contains gas and dust. The number 
density of the spiral population, and the properties of spiral disks 
and the dust therein, are assumed to be either constant in time, 
or to evolve according to evolutionary models taken from the 
literature. 

We take the observed QSO luminosity function as parame- 



terised by Pei (1995) either as it is, or theoretically extrapolated 
towards high redshifts (Haiman & Loeb 1998). 
We can summarise our results as follows: 

1 . At 28 th apparent magnitude, the number density of QSOs 
lensed by spiral galaxies reaches almost 100 per square de- 
gree in the V band, irrespective of whether the observed or 
the theoretically extrapolated QSO luminosity function is 
assumed. The number density of lensed QSOs per square 
degree in the L band is just above 10 for the observed QSO 
luminosity function, and well above 100 for the theoreti- 
cally extrapolated luminosity function. 

2. Neglecting dust evolution reduces the number density of 
lensed QSOs brighter than 28 th magnitude by a factor of 
« 3 in the V band, and by a factor of w 1.5 in the L band. 

3. Neglecting the cosmic evolution of the spiral population 
reduces the number density of lensed QSOs brighter than 
28 th magnitude by a factor of « 5 in the V and the L bands 
alike. 

4. The time-delay distribution peaks at small time delays be- 
low one hour It is bimodal with a secondary maximum 
around one day. The short time delays are due to im- 
ages straddling spiral disks, the long time delays are due to 
image splitting by the halo. Consequently, the secondary 
maximum becomes shallower or higher when the disks be- 
come more or less dominant, respectively. 

5. The fraction of image systems with three or more images 
is of order 30% in the V band and 40% in the L band for 
evolving, maximal disks. It decreases when the disk frac- 
tion is reduced, and increases when the disk evolution is 
neglected. 

6. The redshift distribution of the lensed QSOs peaks around 
z « 3 in the V band, regardless of the QSO luminosity func- 
tion assumed. The QSO redshift distribution in the L band 
also peaks around z = 3 if the observed QSO luminosity 
function is assumed, but the peak broadens and shifts to- 
wards z = 5 — 6 for the theoretically extrapolated QSO lu- 
minosity function. 

These results demonstrate that lensing of high-redshift QSOs 
by spiral galaxies will provide a wealth of information once an 
angular resolution of < 0. 1" can be achieved between the visible 
and the near infrared at flux limits of < 28 — 30 th magnitude. 
The projected NGST can reach such flux limits within at most a 
few minutes of exposure time. The total number of lensed QSOs 
per square degree reflects the number density evolution of spiral 
galaxies out to redshifts of 0.5 — 0.6. The ratio of the lensed 
QSO number counts in the V and L bands is a measure for the 
dust evolution in the spiral disks. The time-delay distribution 
and the fraction of image systems with three or more images are 
sensitive to the disk mass and the scale evolution of the spiral 
disks. 

The disks in spirals are of particular interest here. Spiral 
haloes alone have substantially lower lensing capabilities than 
ellipticals because of their lower velocity dispersions. A sub- 
stantial fraction of the strong-lensing cross sections of spirals is 
contributed by near-edge-on disks. If they are sufficiently mas- 
sive, they give rise to strongly elongated caustic curves with 
naked cusps. Lenses with such caustics produce characteris- 
tic image configurations, namely triple images of comparable 
brightness with very small angular separations straddling the 
disks. The occurrence and abundance of such image configu- 
rations typical for lensing by projected disks will therefore be a 
direct measure for the disk fraction in spirals. 
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